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Abstract. Aims. The spatial distribution of galaxies in the Local Group (LG) is the footprint of its formation 
mechanism and the gravitational interactions among its members and the external massive galaxies or galaxy 
groups. Recently, Pasetto & Chiosi (2007), using a 3D-geometrical description of the spatial distribution of all 
the members of the LG (not only the satellites of the MW and M31) based on present-day data on positions 
and distances, found that all galaxies (MW, M31, their satellites, and even the most distant objects) are confined 
within a slab of about 200 kpc thickness. Examining how external galaxies or groups would gravitationally affect 
(and eventually alter) the planar structure (and its temporal evolution) of the LG, they found that the external 
force field acts parallel to the plane determined by geometry and studied this with the Least Action Principle. 
Methods. In this paper, we have thoroughly investigated the role played by the tidal forces exerted by external 
galaxies or galaxy groups on the LG galaxies (the most distant dwarfs in particular) in shaping their large scale 
distribution. The idea based on the well known effect of tidal interactions, according to which a system of mass- 
points can undergo not only tidal stripping but also tidal compression and thus become flatter. 
Results. Excluding the dwarf galaxies tightly bound to the MW and M31, the same tidal forces can account for 
the planar distribution of the remaining dwarf galaxies. We analytically recover the results of Pasetto & Chiosi 
(2007) and prove that a planar distribution of the LG dwarf galaxies is compatible with the external force field. 
We also highlight the physical cause of this result. 
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1. Introduction 

Over the years, many attempts have been made to find 
the spatial distribution of the LG galaxies. Limiting 
' ourselves to a few pioneering studies and a few con- 
, trib utions in the past decade, Kahn k. Woltier (1959) 
and Ravchaudhurv fc Lynden-Behf i 19891 ) suggested a pla- 
nar distribution based on studies of the LG dynamics; 
lHartwickl ( 2000f ) found a flat ellipsoid with axial ratios 
(a, b, c) — (1.00 0.51,0.19) whi c h is not too di f ferent from 
a plane; finally iKunkell (l979[) . iGrebel et all (|l999t ) and 
iFusi Pecci et al. ( 19951) suggested that the satellite dwarf 
galaxies of the MW and M31 lie on planes 0- 



M ore recently, iKroupa et al. ( 2005 ) and Metz et al 
(|2007l) suggested a planar distribution of the satel- 
lites of the MW, which however could also be ex- 
plained as a conseque nce of the distribution of sub-haloes 



1 Closely related to this problem is the issue of the 
anisotropic distribution of inner sub-haloes with respect to 
large r haloes in relation to the Hol mberg effect dHolmberd 
1969 | ) with diss imilar results, e.g. (|Sales fc Lambasl l2004r . 



Zentner et al. 2005) in the early cosmological stages 
dKang et al.ll2005t). The same problem has been investi- 
gated bv iKoch & Grebel (|2006f ) for the satellites of M31 
with similar conclusions. 

Starting from the basic idea that an off-center hydro- 
dynamical collision occurred some 10 Gyr ago between 
the primordial gas-rich M31 galaxy and the MW, and 
co mpressed the halo gas to form all the LG dwarf galax- 
ies, [Sawa^£ujirngto| (|2005l ) suggested that the new-born 
dwarf galaxies would be located near the orbital plane of 
the MW and M31. They argued that this view is also sus- 
tained by the visual inspection of the 2D sky distribution 



Yang et al.ll2006l ). What matters here (and is still debated) is 



whether disruptions and tidal effects can create the apparent 
polar alignment of the dwarf satellites around the host galaxy 



or, for the particular case of the LG, the position of the dwarf 
galaxies is the consequence of peculiar directions of pre-existing 
cosmological filaments. 
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of the LG members and that a well-defined plane of fi- 
nite thickness is found, within which most of the member 

gal axies are confined . 

iPasetto fc Chiosil (|2007l hereafter PC07) attacked the 
problem from a completely different perspective. In sum- 
mary, adopting known data on positions and distances 
and making use of analytical geometry, they looked for the 
plane that minimizes the distances of all galaxies in the LG 
to it (not only the MW and M31 and their satellites, but 
also the distant dwarfs). The second part of their study 
was to find a dynamical justification for the planar distri- 
bution. To this aim, they applied the Hamilton Method 
(Minimum Action) to investigate the dynamics of the LG 
complex and the action of the gravitational forces exerted 
by external nearby galaxies or groups. They found that 
the planar distribution is fully compatible with the mini- 
mum action and that the external force field is likely com- 
patible with the plane. Such a field pulls the LG galaxies 
along, without altering their planar distribution. Special 
care was taken to evaluate the robustness of the result. 

To so mehow accoun t for t he different resu l ts ob - 
tain ed bv iKroupa et all (120051) . iKoch fe Grebel (|2006t) . 
and Sawa fc Fuiimotol ~ 2005 ). it is worth recalling here 
that the various studies did not use the same galaxy sam- 
pling, start from the same working physical hypotheses 
nor deal with the same dynamical regime. In brief: 

(i) The planes for the MW and M 31 satellites 
(jKroupa et alj 120051: IKoch fe Grebel 120061 respectively) 
are of a very local nature as they are the consequence of 
strong collisional dynamics with the host galaxy (hereafter 
HG). No easy explanation can be found to secure that 
these planes will survive for long time (more than a few 
dynamical time scales) due to the peculiar proper motions 
of the d warfs that determine these planes: see, e.g., the 



ideas injLyndcn- Bell dl983h : iLvnden-Bell fc Lvnden-Bel 



19951) : IPalma " et aiT(|2002l ) applied then bv iMetz et al] 



20081 ). What would be role of distant dwarfs in determin- 



ing the structure of the whole LG is simply not cosnidered. 
Therefore, these planes cannot be extended to the whole 
LG. 



(ii) In lSawa fc Fujimotd (|2005h the solution for a com- 
mon plane is based on the ad hoc initial hypothesis con- 
cerning the origin of the angular momentum. The sample 
of dwarf galaxies used to determine the plane is limited to 
the satellites of the two HGs. Finally, the orbits of these 
dwarfs are constrained to lay on this plane. 

(iii) In PC07, the common plane is chosen by assum- 
ing that it contains the MW and M31 and minimizing 
the distances of all remaining LG galaxies to this plane, 
including also the distant ones. This can be justified by 
considering that the HG satellites are strongly influenced 
by local dynamics (with continuous modification of their 
orbits, including possible captures by the HGs) and that 
if a common planar distribution for all LG galaxies exists, 
this should be brought into evidence by the more external 
galaxies. They are much less likely to be affected by strong 
interactions with one of the HGs and therefore more sen- 
sitive to the influence of external galaxies and/or groups. 



The plane found by PC07 is actually a slab of about 200 
kpc thickness, i.e. it is worth mentioning that the largest 
apo-center of the HG satellites most probably falls inside 
this slab. 

In addition, PC07 have shown that the external force 
field runs parallel to their plane. It is likely that among 
the galaxies of the LG those that feel the external action 
the most are the dwarfs not tightly bound to any HG. In 
other words, this group of dwarfs could act as a tracer of 
the external force field. Our aims is to find a gravitational 
action that is able to induce, on a long time scale, a sort 
of extended slab. Tidal forces are known to engender this 
kind of response. 

The main goal of this study is to highlight the physical 
nature of the results obtained by PC07 that had a rather 
complicated dynamical description requiring a numerical 
approach. We therefore develop here a simpler linear ap- 
proximation that is much easier to handle and yet able to 
provide a physical insight. 

The plan of the paper is as follows. In Section[2j we de- 
fine the tidal force field acting on the whole LG. The tidal 
forces are those developed by external groups of galaxies. 
We lump together the MW and its satellites (M31 and its 
satellites) for which a suitable treatment is required and 
look at the remaining dwarf galaxies. On a long time scale 
the external tidal forces can engender a planar distribution 
of these dwarf galaxies of the LG. In Section[3]we go deeper 
into this issue and check whether the planar distribution 
is a mere coincidence or the consequence of fundamental 
laws of mechanics. The answer is the latter: the planar 
distribution corresponds to a minimum energy and sta- 
ble configuration of the whole system. The dwarf galaxies 
must lie on a plane as a consequence of the long time scale 
influence of the tidal forces exerted by massive galaxies or 
galaxy groups external to the LG. Furthermore, the plane 
found by PC07 and the minimum energy plane are coin- 
cident and the situation is stable. Finally, in Section [4] we 
summarize the results and present some general consider- 
ations. 



2. LG structure and tidal forces 

Looking at the composition of the LG, three main com- 
ponents can be identified: two massive galaxies (MW and 
M31), their respective groups of bounded satellites, a large 
number of distant dwarfs galaxies loosely interacting or 
even uncoupled to the dominant galaxies. As far as the 
gravitational interaction is concerned, there are two ques- 
tions we are interested in addressing: 

(1) How does the external force field change with time? 

(2) Have the dwarf galaxies that are not members of 
the MW or the M31 family been affected by the tidal 
interaction with the external force field? 

First, to get a rough estimate of the influence of the 
tidal forces acting on the LG, we must develop an accu- 
rate geometrical description of the LG during its temporal 
evolution. Although the tidal forces are weak, they can 
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Table 1. External galaxy groups gravitationally influenc- 
ing the LG dynamics. Each group is indicated by the dom- 
inant galaxy. The members of each group are the same as 
in PC07 to whom the reader should refer for all details 
(their Section 4 and Table 2). To this list the MW and 
M31 are added. The radial velocities are quoted relative 
to the center of the MW. No correction for the motion 
of the Sun toward the Local Standard of Rest is applied 
because it falls below the accuracy adopted in this study. 
The uncertainties on the distances and radial velocities 
are omitted. 



Group 


I 


b 


d 


M 


V, 


Name 


O 


O 


Mpc 


1O 12 M 


km s -1 


IC 342 


138.2 


10.6 


3.3 


12.6 


171.0 


Maffei 


136.4 


-0.4 


3.5 


6.3 


152.0 


M 81 


142.1 


40.9 


3.7 


1.6 


130.0 


Cen A 


309.5 


19.4 


2.7 


4.7 


371.0 


Sculptor 


105.8 


85.8 


3.2 


6.3 


229.0 


M 83 


312.1 


25.9 


4.5 


0.8 


249.0 


Andromeda 


121.21 


-21.60 


0.76 


3.16 


-123.00 


Milky Way 








0.00 


2.20 


0.00 



produce a cumulative effect of compression that could ex- 
plain the planar distribution by acting over a very long 
time. However, if this action is not always pointing in the 
same direction, the net effect can be small, and the tidal 
effect can no longer be the physical mechanism compress- 
ing the dwarf galaxy distribution. To follow the direction 
of a tidal fi eld, we apply t he us ual formalism of the tidal 
tensor (e.g. Misner et al. ( 19731 ). Chap. 1) to all galaxies 
of the the LG and try to understand its global behavior 
under the action of the external force field. 

Frame of reference. The actual mutual interaction 
between the MW and M31 suggests that they can be con- 
sidered as a privileged system whose center of mass (CM) 
can be assumed as the origin of a reference frame (pre- 
sented in Fig. [T] and described in more detail below). The 
positions and motions of any other dwarf of the LG that is 
not a member of the MW and M31 families can be given 
relative to this CM. This is reminisce nt of the g e ometr y 
of the restricted 3 body problem (e.g. Szebehelv ( 19671 )) 
but here the CM will be followed in its time evolution. In 
particular the stationary action principle has been applied 
in PC07 to produce a possible solution for the motion of 
the external groups acting gravitationally on the LG as 
well as MW and M31 (their table 4). From this solution 
we can infer the spatial evolution in the time, t, of the 



entire LG-barycentric system x 



CMlg 



„CM LG 



(i), which 



is nearly coincident with the center of mass of the MW + 
M31 system. 

Here, we start by considering the general expression 
for the external tidal tensor due to any potential $ ex- 
pressed in the reference system So of Figure Q] to express 
the external tidal tensor acting on the complex MW+M31 



( X CM LG ) = 



GAL 




Fig. 1. Sketch of the geometrical framework we have 
adopted. First we define the inertial reference frame, al- 
ways named So, with axis (So, X\, x%, x%). It is centered 
on the CM of the external galaxy groups listed in Table [TJ 
We then introduce two auxiliary reference frames: the first 
centered on the barycenter of the MW and M31, which is 
aligned with the principal axes of the eigen-system pro- 
vided by Eqn. fl}. This system is called Si and it has 
axes (£1,^27 £3) with £ 3 pointing roughly in the direction 
orthogonal to the plane of the dwarf galaxies (red dot- 
ted rectangle). The second is called S2 and it has axes 
(Vi j V2, with 773 pointing in the same direction of £3 and 
j]2 pointing in the direction of M31; see the text for details. 
The current position vector of a generic dwarf galaxy in 
the system Si and S2 is called £ and rj respectively (e.g. 
we show £ in the figure). It is always oriented toward the 
generic dwarf galaxy. 



/3 



00 



.CMlg 
j 



X- 



V 



CMlg 



(1) 



where G is the gravitational constant, x s are the coor- 
dinates of the external galaxy groups of Table [T] of mass 
M g , 5 is the Dirac's function, ||.|| is the standard norm. 
We evaluate the tidal tensor at the barycenter of the LG, 



g#MW,M31 xf MljG - xf 



x cm lg _ jjUMlg ^ ag (i emiec i above but where the time 
dependence has been omitted. 

To get an idea of the effects we are looking for, let 
us make the following preliminary considerations. Let us 
assume that the most distant dwarf tat define the plane 
7r discovered by PC07 owe their distribution to the ini- 
tial conditions determined by the tidal tensor of exter- 
nal objects. Then we expect the total tidal tensor (TTT), 
X|j 0t (x), defined by Eqn. |T]) including in the sum also 
the MW and M31, to give, if reduced to its normal form, 
the most negative eigenvalue. Between the three eigenvec- 
tors of TTT, the one corresponding to this most negative 
eigenvalue will indicate the direction of the tidal compres - 
sion (se e classical text-books such as iMisner et aL ( 1973 ) 
Chap 1, Binnev fc Tremaine ( 19871) Chap 7. or som e appli- 
cations as in IRavchaudhurv &: Lynden-Bell (|l989fn . Thus 
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Table 2. This table shows the temporal evolution of the eigenvalues and eigenvectors as a function of the look-back 
time in Gyr (column 1). The three eigenvalues are in columns (2) to (4) and the corresponding eigenvectors are in 
columns (5) through (13). The headers are self explanatory. We have highlighted the most negative values (column 
4) and the corresponding eigenvectors (column 11 through 13) for all values of the look-back time at which the most 
negative eigenvalue keeps its sign. 



lib 


M 




A3 


n x 


fly 


Tlz 


„(A 2 ) 
«x 


Thy 


„(>*) 
n z 


„(A 3 j 

n x 


Thy 


Tl z 


13.3 


0.0125 


-0.0049 


-0.0076 


0.310 


0.507 


0.804 


0.934 


-0.008 


-0.356 


0.174 


-0.862 


0.476 


12.3 


0.0058 


0.0001 


-0.0060 


-0.266 


0.411 


0.872 


0.944 


0.292 


0.151 


0.192 


-0.864 


0.465 


11.0 


0.0031 


0.0008 


-0.0038 


0.941 


0.052 


-0.334 


0.265 


0.498 


0.826 


0.209 


-0.866 


0.455 


9.5 


0.0030 


-0.0002 


-0.0028 


0.971 


0.235 


-0.043 


-0.066 


0.436 


0.897 


0.230 


-0.868 


0.439 


7.9 


0.0030 


-0.0006 


-0.0024 


0.954 


0.297 


0.026 


-0.146 


0.390 


0.909 


0.259 


-0.871 


0.416 


6.2 


0.0030 


-0.0007 


-0.0023 


0.939 


0.341 


0.045 


-0.171 


0.346 


0.922 


0.299 


-0.873 


0.384 


4.6 


0.0033 


-0.0008 


-0.0024 


0.919 


0.391 


0.050 


-0.175 


0.292 


0.940 


0.353 


-0.873 


0.336 


3.0 


0.0037 


-0.0009 


-0.0028 


0.887 


0.458 


0.059 


-0.171 


0.208 


0.963 


0.429 


-0.864 


0.263 


1.4 


0.0046 


-0.0010 


-0.0036 


-0.829 


-0.551 


-0.097 


-0.158 


0.065 


0.985 


0.536 


-0.832 


0.141 


0.0 


0.0065 


-0.0017 


-0.0047 


-0.738 


-0.657 


-0.151 


-0.064 


-0.154 


0.986 


0.671 


-0.738 


-0.071 



if we find this behavior also in our TTT evaluated at the 
position of a dwarf galaxy far away from the MW or M31 
this could hint that the compression docs indeed occur. 
For example a simple case would be that a dwarf presently 
belonging to the plane ir was also formed by some mech- 
anism in the plane n or close to it. In such a case we can 
simply make the hypothesis that the tidal force acting on 
this dwarf should be similar to the tidal force that acted on 
the plane in the past, say 9 Gyr ago. We can estimate these 
eigenvalues and their eigenvector directions by combining 
the equations of the plane ir (see Eqn. ((2])), obtained by 
looking at the current dwarf galaxy distribution, and Eqn. 
|T]) which can also be evaluated backwards in time. If we 
find compatible values between different points at different 
epochs, then we can claim that the effect we are searching 
for could effectively have acted. 

If the plane ir is a slab with a diameter of 4 Mpc 
and 200 kpc thick, we can evaluate the TTT at any 
point on this plane P G n, say 2 Mpc away from the 
barycenter of the LG, 7]j 0t (x p ). The resulting eigen- 
value of this tidal tensor, e.g. 9 Gyr ago, is Ag Gyr = 
{0.024, -0.016, -0.007}. The same evaluation can then be 
repeated for the barycenter position Ty ot (x CMlg ) , obtain- 
ing Ag Gyr G = {0.022, -0.014, -0.007}. This result strongly 
suggests that back in the past the force determining the 
subsequent orbital evolution of a generic dwarf had a com- 
ponent squeezing the motion toward the plane. Proceeding 
in this way we can prove the compatibility of the eigen- 
values of the TTT for every position on the plane 7r, i.e. 
j^tot ( x CM LG j ^ jrtot ( X P) VP g 7T. This clearly allow us 
to explore the possibility that the plane ir of PC07 is the 
consequence of the tidal forces acting on the LG during a 
large fraction of the Hubble time, i.e. we want to extend 
this estimation not only to the present time t = to but 
also to the tim<H t < t . 



2 Of course a more correct computation could have been per- 
formed by knowing distribution of the dwarf galaxies in the 
past, but unfortunately we cannot track back the past orbits 
of the dwarf galaxies belonging to the plane n today; therefore 



To proceed further, we need to search the eigenvec- 
tors associated with the tidal tensor with the most nega- 
tive eigenvalues. In the limits of our approximation, they 
should keep a direction with respect to an inertial refer- 
ence frame not too far from the normal to the geometri- 
cally plane h (it) for most of the Hubble time, which can 
nowadays be inferred by simple inspection of the dwarf 
galaxies' distribution in the LG. To prove this, we solve 
the eigen-system for the tidal tensor of Eqn. ([1]) as a func- 
tion of time. The solutions are presented in Table O The 
results we are interested in are limited in time to a range 
where monotonicity of the trend of the eigenvalues can be 
exploited in order to reveal an integrated cumulative ef- 
fect of compression or expansion. We find that our range 
of interest has to span the last 9 Gyr, imposing a lower 
limit to our analysis of t = t m { = 9 Gyr. Before this £ in f the 
configuration of the eigenvalues is slightly different. From 
Table [2] we see that the time evolution of the eigenvalues 
shows a phase with one compressive direction compared to 
two positive expansion directions, and before that, as well 
as after t = i; n f, we see two negative directions compared 
with one positive. For simplicity we will not be treating 
analytically in the following sections these switchings be- 
tween the configurations. We are only interested in the last 
most dominant time evolution of the monotonic behavior 
of the eigenvalues. As we will exclude the first three rows 
of Table [2] (the primordial evolution prior to ii n f) from our 
analysis, from now on we will uniquely identify (unless 
otherwise specified) with Ai, A2, A3 the positive eigen- 
value, the second negative eigenvalue and the most nega- 
tive eigenvalue respectively for the eigenvectors of the ex- 
ternal tidal tensor defined in equation {1} followed in their 
time t evolution exclusively for t e ]t- m {, to}. The most neg- 
ative eigenvalue (column 4) and its evolution during the 
past 9 Gyr and projections of the associated eigenvector 
onto the axes of the inertial system are highlighted in ital- 



Eqn. ((2]) for the plane 7r cannot be directly determined in its 
time evolution. 
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0.015 




Fig. 2. Evolution of the eigenvalues #(T;? xt ) of Tabledas 
a function of the look-back time in Gyr (column 1). The 
red line with dots is for the first positive eigenvalue Ai, the 
blue line with squares is the second eigenvalue A2 negative 
for tib < 9.5Gys, and the black with triangles is for the 
third (most negative) eigenvalue A3. 




2 4 6 8 10 12 14 



Fig. 3. Here we tracked back the influence of each indi- 
vidual component in the sum of Eqn. (fTJ) for the most 
negative eigenvalue, A3' . The curves are normalized to the 
black line, A3, of Fig. [2] (see text for details). The ma- 
jor contribution overall is dominated by the IC342 group, 
followed by Cen A, Maffei and Sculptor groups. The con- 
tribution of M83, even if nowadays marginal, could have 
played a role in the past in the imprint of the initial con- 
dition of the LG dwarf galaxies. 



ics. The time variations of the three eigenvalues are shown 
in Fig. H 

It is soon evident from the data displayed in Fig.[5]that 
one of the negative eigenvalues dominates. Therefore, we 
are in a situation in which the external field acting on 
the LG may engender a planar distribution of the dwarf 
galaxies. This is an interesting result because 
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r, we can analyti- 
cally obtain the same results of PC07 for the behavior 
of the external field. The external force field turns out 
to be compatible with a flat spatial distribution of the 
dwarf galaxies that remains constant for a large frac- 
tion of the Hubble time. Moreover, we can argue that 
the external force field started to flatten the spatial 
distribution of the dwarf galaxies already prior to ij n f . 
If we investigate in more detail the relevance of the dif- 
ferent groups on this flattening effect, we can plot in 
Fig.[3]the normalized trend of the most negative eigen- 
value of Eqn.fT]) split into its components. Here we see 
the most negative eigenvalues of the sum in Eqn. (fTJ) 
normalized to the overall sum (hence the black line of 
Fig. [2] is here the unitary constant upper bound of the 
figure). As we can infer from the figure, the influence 
of IC342 group has always been the most significant, 
followed by the effects of Maffei, Sculptor and Cen A 
groups. This is expected from their masses and posi- 
tions listed in Table [TJ The time evolution shown in 
this Figure confirms their relative importance in the 
compressing effect on the LG for its temporal evolu- 
tion. Slightly less important is the contribution of the 
M83 group that was nevertheless as important as that 
of the M81 group 11 Gyr ago. 

2. We must clarify once and for all that the present result 
does not prove that the spatial distribution of dwarf 
galaxies in the LG has to be flat, but only that the 
external force field is compatible with such a flat dis- 
tribution. 

3. The fact that the external force field acting only on the 
two main HGs (MW and M31) of the LG is compatible 
with a flat spatial distribution of the distant dwarfs 
does not tell us anything about the distribution of the 
nearby dwarf satellites around their HG (both MW 
and M31). 

4. Furthermore, the planar distribution shown by the geo- 
metrical analysis made by PC07 cannot a priori be re- 
lated to the planar distribution suggested by the com- 
pression effect described by the tidal tensor affecting 
the LG. The subject of the following analysis is to ex- 
plain the coincidence claimed by PC07 between the 
geometrical and dynamical planes. 

Now we seek to prove that the following two issues are 
tightly related: (i) There exists a plane in the spatial dis- 
tribution of the LG dwarf galaxies that is expected from 
the external tidal force acting on the LG. (ii) In the con- 
text of the linear approximation that we have adopted (see 
below) , the tidal force field compatible with a flat distribu- 
tion has an orientation whose normal vector is compatible 
with the normal vector of the planar distribution 7r found 
by PC07. 

The equation for the plane tt can be rewritten here in 
the reference frame So as 

0.64a; 1 - 0.61x 2 - 0.45x 3 = (2) 



6 



S. Pasetto & C. Chiosi: Tidal effects on the Spatial Structure of the Local Group 




30 E 60 E 90 E 120 E 150 E 180 E 



definition of the direction of h (tt) , the normal to the ge- 
ometrical plane and the direction of the vector associated 
with the most negative eigenvalue are compatible at the 
present time at the 3a-level of confidence. Moreover, the 
eigenvector of the most negative eigenvalue A3 of the tidal 
tensor moves, maintaining a direction not so far form the 
direction that we can nowadays deduce from the observa- 
tion for n(7t). This result, which was already present in 
PC07, is recovered here in a semi- analytical treatment of 
the whole problem. The key question to be answered now 
is whether the coincidence is "causal or casual" and what 
the meaning of all this is. 



Fig. 4. The Hammer projection of the sky position of the 
normal to the plane found by PC07 (from purely geomet- 
rical arguments). The color code from light red to dark 
brown corresponds to the 3<7, 2er, and la uncertainties as 
estimated by PC07 from the principal component anal- 
ysis and projected onto the sky. The shape of the three 
shaded areas becomes more oblate at increasing angular 
direction b. The classical Tissot's circles, perfect circles of 
angular radius of 7 deg have been plotted to help visual- 
ize the angular distances. We have used the same symbols 
and time intervals as in Fig. [2] and Table [2] The black tri- 
angles nearly overlap the shaded area, showing that the 
direction of the eigenvector with the most negative eigen- 
value lies close to the direction of the normal to the PC07 
plane, which is of course constant in the figure because it 
is derived from the observational data available today. For 
comparison we also show the evolution of the other two 
orthogonal eigenvectors indicated by the same symbols as 
in Figured 



with a direction (l,b) = (45,27) for n(7r), the normal 
vector to the plane. Furthermore, the directions of the 
three basis vectors es in which the eigenvalues of Table [2] 
are expressed have been chosen as collinear with those of 
the reference frame adopted by PC07 (the approximation 
was made such that the Sun is placed at the center of MW, 
the error of « 8.5 kpc is negligible in the present context). 
Therefore, we can assume that only at the present time 
does the reference system in which the above equation 
of the plane is written have orthonormal vectors parallel 
to the basis vector of reference system adopted at t = to- 
It follows that the Hammer projection from the center 
of mass of the LG shows the eigenvector directions as a 
function of time and the normal vector to the geometrical 
plane of PC07 are as in Fig. [2l 

Interestingly, it is the eigenvector relative to the most 
negative eigenvalue that lies closer to the direction of the 
normal of the plane of the dwarf galaxies. 

Considering that every point in this map has an error 
radius of about ±9°, inherited from the Minimum Action 
analysis, and considering the uncertainty in the angular 



3. Causal or casual? 

The tidal tensor can be derived from the Taylor expansion 
of the force field (see next section). This implies that the 
object we want to investigate (a dwarf galaxy) is under 
the effect of a smoothly varying potential in the course of 
evolution. 

The large scale description of the gravitational interac- 
tion adopted here and by PC07 does not work at the dis- 
tance scales of the closer HG-dwarf satellite interactions 
where satellite dwarf galaxies suffer kick-off, are continu- 
ally absorbed into the halo of the HG (MW and M31 in 
our case) and the collisionless description is not correct. A 
dwarf galaxy of these closer samples, undergoing the much 
more intense direct interaction with the HG, does not sig- 
nificantly respond to the weaker external field acting on it. 
On the other hand, the collisionless description is suited 
to deal with the effects of distant galaxy groups (see the 
list in Table [T]) and hence the influence of the external 
field on the other dwarf galaxies fa r away from a HG (se e 
e.g. iRavchaudhurv fc Lynden-Belll dl989n : IPeebled (|l990h : 
iDunn k Laflammel (|l993h : IPeebled (|l994h ). 

The frame of reference. In order to formalize what 
is in the previous section we present the following two 
hypotheses, one on the directions of the eigenvectors and 
the other on their values: 

1. We will not consider the evolutionary stages of the 
Universe older than tvo > 9 Gyr (look-back time) to 
work with a monotonic behaviour of a single negative 
eigenvalue. 

2. We assume that the eigenvector corresponding to the 
most negative eigenvalue always points in the same 
direction with respect to the inertial reference frame 
in the non-comoving coordinate system. 

First we define Sq as the inertial reference frame cen- 
tered on the CM of the external galaxy groups listed in 
Table [1] We introduce also a second reference frame Si , 
non-inertial and comoving with the CM of the LG, ap- 
proximately coincident with the CM of the pair M31 and 
MW. Moreover, we assume that this reference frame has 
its principal axis always collinear with the principal axis of 
the external tidal tensor, i.e. this frame is an eigen-system 
of the external tidal tensor. XM31 is the coordinate vector 



S. Pasetto & C. Chiosi: Tidal effects on the Spatial Structure of the Local Group 



7 



of M31 in Sq, xmw is the coordinate vector of MW in So, 
x is the vector of a generic dwarf galaxy in So and R is 
the position of the M31-MW barycenter, i.e. the origin of 
Si in Sq. The generic position vector in the frame Si is £. 
It will be used to indicate the generic position of a dwarf 
galaxy. Both reference frames are orthogonal (see Fig. [T|) . 

Equations of motion. We only point out here that 
the generic position vectors in the two reference frames 
are now written as x = x™ + 0£ = R + 0£, where 
O 6 SO (3) is the generic rotation matrix with elements 
Oiy Similarly, for the time derivative of x we can write 

x = ^lg + O ( £ + ^ x £ ) or, in more compact notation 



where ft is the angular velocity and we put V = x 



CM 
LG ■ 

The time derivative of the velocity equation yields the 
equation of motion 



A + o(|+2nx£ + nx£ + nx(nx£ 



where, in the second equation the vector A has been 
used to simplify the notation for the acceleration x^q oi 
the barycenter. Moreover, to better underline the physical 
meaning of the different terms we recollect the previous 
equation as 



together with two other similar relations for the poten- 
tials of MW and M31 that are not given here for the sake 
of brevity, with (., .} expressing the standard inner prod- 
uct. Inserting these three relations in Eqn. Q, after some 
simplifications we obtain the acceleration x — A in linear 
approximation: 

x - A ~ (T ext (R) + T MW (R) + T M3 i (R),x-R). 

By introducing the TTT as defined in the preceding sec- 
tion with 

T (R) = T cxt (R) + Tmw (R) + 2m31 {R) 

and using x — R = 0£, we can simplify Eqn. (J3J) that 
becomes 



£ = O t TO£ - 2SI x £ - SI x £ - SI x (ft x £) , (6) 

where we have dropped the explicit dependence on the po- 
sition of the tidal tensor T. This tensor is always evaluated 
at the barycenter T — T (R) if not specified otherwise. 



3.1. Energy states of equilibrium 

Our aim now is to understand whether the equation of 
motion ^ can lead to stable equilibrium configurations. 
There are several techniques for attacking this problem 
based on the integration over a time interval of th e force 
or the impulse, see e.g. ( Binnev fc Tremain^ 19871 Chap 



£ = O (x — A) — 2SI x £ — SI x £ — SI x (SI x £) . (3) 7), th e elegant Lagrangian treatment of iGnedin et al 



From this last equation it is quickly evident that the 
acceleration £, suffered by a dwarf galaxy in the non- 
inertial reference frame, is the sum of different terms due 
to the motion of the frames Si and So and the fictional 
forces that appear thanks to the non-inertial nature of 
Si: — SI x (SI x £) is the centrifugal effect, — 2SI x £ the 
Coriolis effect proportional to the velocity of the dwarf, 
and — SI x £ is the effect caused by the non-constant ro- 
tation rate of Si. As usual we assume that the rotation 
matrix linking the different orthonormal reference frames 
is O : OO t = I where the T stands for transpose and J 
is the identity matrix. 

From Eqn. ([3]) it is easily evident how to evaluate the 
term which gives the acceleration of a generic dwarf, x— A. 
Now the acceleration of a dwarf galaxy in So is due to 
three contributions: the gradient in the external potential 
V$ ox t, the gradient in the MW potential V$mw, and the 
gradient in the M31 potential, V$M3i, i-e. 



[1999) , or the soph isticated analysis in the action space 
of lWeinberd (ll994al lbT). 

To proceed further we look for the energy equilibrium 
configurations and their stability of a system governed by 
the equations of motion ([6]). Here, we take advantage of 
the fact that we can follow the evolution of the angular 
velocity of the non-inertial reference frame by looking at 
the motion on the sky of the eigenvector of the tidal tensor 
as already made in previous sections for the quadrupole, 
see Fig. [3] 

It can be demonstrated that the Lagrangian leading to 
Eqn. up to a total derivative, can be written as 



L 



m 
~2 



£,ttx£ 



m 
~ ~2 



|n x £|| 2 -m(A,0£}-W 



(see for instance Landau fc Lifshit3 1969h . Therefore, tak 



X = -V x ($ e 



MW 



M31 



(4) 



ing the linear momentum p = ^ 
can evaluate the energy E = (p, £ 



E 



Taylor expanding <j> cxt at the second order around the 
comoving barycenter, and with the usual definition of the 
Tidal Tensor given in equation (TTJ), we can easily write 

$ cxt (x) ~ $ oxt (R) + (Vi^cxt (R) ,x — R) — 

1 (5) 
- T ext (R) (x — R) , x — R) 



m£ + m (SI x £) we 
-I as 



-y \\Slx$f +m(A,Ot) + W. 

Taking the derivative with respect to the positions, we 
obtain 



dE 



mSl x (SI x £) - mO T TO£. 



(7) 



The equilibrium energy is then given by the solution of 
the equation 

SI x (SI x £) - O t TO£ = 0. (8) 
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In the reference frame Si, by definition the external tidal 
tensor is always diagonal, i.e. T T cxt (R) O is greatly sim- 
plified but, in contrast, the first term of Eqn. © has a 
complicated structure and the same occurs for the other 
two terms composing T. This way of proceeding does not 
allow for significant simplifications but a better insight 
can be gained by moving to a new reference frame, S2, 
with the same origin as the previous reference system Si , 
in which the matrices representing the tidal tensors Tm3i 
and Tmw are both in diagonal form, but where T ex t is 
not. From the analysis of the orbit evolution expected in 
a statistical interpretation of the Minimum Action, PC07 
showed that MW and M31 are roughly coplanar to the ex- 
ternal force field, at least during the last 9 Gyr. Therefore, 
in our simple model we can assume that the angular ve- 
locity vector fijg with which we describe the rotation of 
the system of reference Si , which we remember here again 
for sake of clarity has its third axis £3 parallel to the nor- 
mal n (tt) , has to be parallel to the angular velocity vector 
in this new system of reference (O, e m , e m , e^), namely 
0?) 3 (S2), that we assume rotating around its third axis 
773, i.e. n v3 (S2) 1 1 (Si) over the last 9 Gyr. This means 
that we do not allow the orbital plane of M31 and MW to 
tilt with respect to the plane t\ 3 = it orthogonal to the 
eigenvector corresponding to the most negative eigenvalue 
A3 (t) in the course of time evolution. Even though this as- 
sumption is automatically fulfilled in the analysis below, 
before looking at the evolution of the MW and M31 in the 
sky centered on the LG barycenter, we have checked that 
the angular distance between the position vectors of MW 
or M31 and h (t\ 3 ) is 90°±3°, thus indirectly showing that 
the propagation of numerical errors in the calculation of 
the temporal evolution of the eigenvectors is small. 

On the basis of these considerations, we can suppose 
that 3 N = N (t) : £ = Nrj, i.e. there exists a linear oper- 
ator represented by a rotational matrix TV £ 5*0(3) such 
that the generic position vector of a dwarf galaxy in Si, £, 
can be written as a function of the generic position vector 
in S2, rj, and the tidal tensor of the external potential can 
be written as: 

' Ai \ / Nji N12 \ (m 
T T ext ON V = I A 2 \ [ N21 N22 ?72 

A 3 / V 1 / V % 

7Vii?7iAi + Ni 2 r/ 2 Xi 
N 2 ir]i\2 + /V22772A2 
773A3 

where we have adopted a rotation matrix spinning about 
fi^3 . As we can see, no explicit dependence on the ro- 
tational coefficients of the matrix O is necessary in this 
reference system 5*2 that we have chosen but the reader 
should keep in mind that O (t) ^ N (t) Vi and only in this 
new reference system £2 have we been able to simplify the 
matrix as above. The general form of the rotational ma- 
trix could be casted as a combination of trigonometric 
functions, but this would be superfluous here. The only 
important thing to note is that even if the assumption of 
collinearity between the axis of rotation of Si and S2 has 



been justified, nothing can said about the moduli of their 
angular velocities. We cannot assume that the rotation 
of the pair MW and M31 or, equivalently, the reference 
frame tightened to this rotation spins with the same ro- 
tational velocity of the reference frame attached to the 
external potential! This could lead to wrong or paradox- 
ical results that need to be avoided. In other words we 



cannot assume ||S~2 r?3 (S2 



|f2j 3 (Si)||, which is wrong, 



but simply adopt the condition on the directions 



fl m (S2 



"fe (Si 



\\n 713 (s 2 )\\ \\n i3 (Si)\\ 

For example we can explicitly write the matrix N as 



Nu N12 
N21 N 22 
1 



cos 7 —sin 7 
sin 7 cos 7 
1 



where 7 = 7 (t) is the angle between the two systems. It 
will change as a function of time due to the time depen- 
dence of the angular velocity of the two frames 



7 = 7 



Slg (t) ,Sl% (t) 



This also means that the angular velocity is not related to 
the angular momentum in a simply way. It will indeed be 
the result of the combined action of the centrifugal force 
due to the rotation of the frame tightened to the external 
force field, i.e. T oxt , and the centrifugal force due to the 
rotation of the frame tightened to the motion of M31 and 
MW via Tmw and Tm31- We have a double centrifugal 
effect: one caused by the external field and the other by 
the MW and M31 that acts with different characteristic 
angular velocity (typically time-dependent angular veloc- 
ities). The intensity of the force due to the centrifugal 
component can be written as 



iivn(Si)ir 
= iifi(Si)ir 

= f 110(50 



||2Vr,|| + ||0(S 2 

NI + l|o(s 2 )|| : 

2 + l|0(5 1 )f 



1*7 1 



WW 



\W\ 



where ||t7|| is the distance of the dwarf galaxy from the 
moving barycenter of S2 • Moreover in this frame of refer- 
ence S 2 we have that T T MW ONr) and T T M3 iONr) 
are diagonal . Ther efore, we immediately have (see e.g. 
Misner et all (|l973l ). Chap 1) 




O t T mw ON V = 



where for simplicity we have assumed that e Tj2 points from 
the barycenter of the system toward MW, e V3 is parallel to 
fi?2 (^1)1 and with the third axis oriented in such a way 
to form a left-handed reference frame. In the same way 
we may write the analogous for T T M 3iONr). Finally, 
the tidal term in the condition for the energy equilibrium, 
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Eqn. 



is given by 



O t TO = 



with 

ail = NnXi 




GMmw 



GMm31 



n^Mwir iiT?M3iir 

i\r \ i 2GA/mw i 2GAfM3i 

0(22 — JV22A2 + j- p- + it- rrr 

||»JMW|| ll»7M3i II 

(Tqq — Aq GMmw GM M 3i 

33 3 ||»JMW|| 3 ||rjM3i|| 3 

which provides a simple description of the tidal effects. 
The above relation can be further simplified by recalling 
that, from the definition of barycenter in the reference 
frame in use, we can write 



MW 
V2 



% M31 M M3 1 



Mi 



MW 



(9) 



Inserting this expression into Eqn. ([5]) 
the S 2 frame, and using the relation 

|3 I MW 
72 



written 

I^MwII 



. 0+(% MW j- + = \^\ i.e. I^Mwir = KT, af " 
ter tedious algebraic simplifications, we get for the three 
components of the tidal energy 



O'TOn 



V 



N 2 i\ 2 T]i + 
T]3 



A, 



iViiAi 
-^22 A2 

G«w+<3i) 



^ 31 h M3i r 

2G(Mmw + Mm3i) 



/ 



In the same way we can derive the term due to the 
apparent force. Switching to the revolving system S 2 we 
apply another rotation expressed by Nflx (NflxNrj) = 
N (fix (ft X T))) where ft = fl (ft (Si) , ft (S 2 )) is the an- 
gular velocity of the system S 2 as seen from Sq. The re- 
sulting centrifugal term is 



ft x (fl x 77) 



(N 2 im 



N 22 m) n§ 







The associated equilibrium energy state is given by the 
solution of the systems 



-N 12 n 2 (Ai 
-27772 - (iV 2 i?/i 
(7 - h)V3 = 



+ Vi (7 - NnAi 
^22%) (A 2 



= 



where 



G M 



7 



MW " 



M 3 , 



,„M31| 
M31 |'/2 



> OVt. 





(10) 

(11) 



The system (|10[) is t/ie result we are looking for. It is 
evident from this system of equations that the plane 773 = 
is the equilibrium plane for the dynamical evolution of the 
gravitational system. 

The above result is fully adequate for our purposes 
because the total potential is separated into the radial 



and vertical components. This is a standard consequence 
of the linear approximation obtained by truncating the 
Taylor expansion of the potential at the second order. 
Equivalently one could use a generating function satisfying 
the Stakel theorem for separability in a Hamilton- Jacobi 
equat ion for the above system (e.g Boccaletti fc Pucaccol 
(|l998h ). 

As far as the temporal evolution is concerned (t G 
]£infi*o])j we have 773 (t) — which, translated into our 
spatial resolution, simply means a physical spatial reso- 
lution 1 773 1 < 100 kpc. Furthermore, one can never have 
(7 — A3) =0 because this would imply that 3t : A3 (i) = 
7 (t), whereas according to definition (TlT]) we have 7 (t) > 
Vi, and finally A 3 < V*n> < 9 Gyr as shown by Fig. O 
These last two conditions are clearly inconsistent leading 
to a contradiction that concludes our proof as required. 

Therefore, the major conclusion of this 
demonstration is that the only possible solution is 
the following one: the statistical minimization of 
PC07 is compatible with a planar distribution of 
the dwarfs and it is not a mere coincidence. This 
result completes the missing interpretation of the result 
already included in PC07. 

3.2. Stability of the equilibrium configuration 

Finally, we examine the stability of the equilibrium plane 
that we have found by solving Eqn. (jlOjl . We take the 
energy of Eqn. written for the system S 2 



dE 
drj 



mN [fl x (fl x 77)] - mO T TONr] 



and calculate the derivative 



d 2 E 



dndn 

where £ijk is the Levi-Civita symbol, summation over 
the repeated index is assumed and we define [rW]y = 
[O t TON] .. . Indicating with S = ON the composition 
of the two rotation matrices O from So to Si and TV from 
Si to 5*2, after some algebraic simplifications we obtain a 
more compact form 



d 2 E 



= -mN 



{® + S T Ts} 



dr/dr/ 

where we have defined another matrix 



(12) 



inir&i 



-Qifi 2 



-ftift 3 
-n 2 n 3 



This matrix has no inverse, determinant |0| =0, and trace 
Tr (6) = 2||fi|| 2 . 

At this stage the usual procedure would be to solve for 
the eigen-system (fT2|) . However, we can avoid this com- 
plication by noticing that in the linear approximation the 
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vertical potential decouples from the radial one. In other 
words, applying the full procedure we would obtain three 
eigenvalues {a\ (t) , at2 {t) , 7 — A 3 (t)}, where oq (t) and 
0.2 {t) are two complicated functions of the time, whereas 
the last eigenvalue has to be exactly the most negative 
eigenvalue of the external tidal field for which we have al- 
ready calculated the time dependence as shown in Fig. [2] 
As long as the eigenvalue remains negative, the relation 
7 (t) - A3 (i) > 7 (t) > A 3 (t) holds because 7 > by 
definition (equation (|10|)) whereas A3 < as required in 
Fig. [21 Therefore the plane tt is stable. 

3.3. Equations of motion and force balance 

Finally, we can get a much deeper insight for the physical 
reasons of the existence of the plane tt by analyzing the 
equation of motion in 6*2. If the plane tt is a stable config- 
uration of the spatial distribution of dwarf galaxies, it is 
important to isolate the force acting on it. The equations 
of motion in £2 are, from e.g. Eqn. J6j, 

77 = N T ((Text (R) + T MW (R) + T M3 i (i?) , Nq)) - 
2tt x 77 — tt x n tt x (tt x n) 

where for the different elements of the tidal tensor we can 
now write 



Ar T T M3 iiVr7 = 



ll*7M3l|| d 





2GMM3HI2 131 



V 







GM M3 iyf 31 

iimmlP 



an analogous equation for MW, and 



N T T oxt N = 



Tn T 12 
T\2 T22 
A 3 



The equilibrium in the tt plane is given by 



V3 



WVM31 II 

-TfrS% 
- {^22 - 

_GMm31 
ll»7M31 f 



ll»7Mw|h 



y GA/m31 GMmw 

T\2T\2 
T12V1 



Vi- 



2n 3 ri2 + n 3m = 

2GA/M3 



■ (I3V1 = 



2 GMmw \ irt 



IWmw II 



+ A 3 = 



which sheds light on what is happening in realitj{j| 

Along the direction orthogonal to the plane tt, three 
forces are present. One is due to the MW, — F MM ^f-i , an- 

other due to M31, - f? MM3 h , and the third one due to the 
external field +A 3 < 0. They all sum together to push any 
dwarf which tends to escape from the position of equilib- 
rium to return back in the plane. Therefore they tend to 



3 We did not exploit here the Eqn. <(9j previously necessary 
for the energy analysis, in favor of a clearer and easier physical 
interpretation of the terms in the equation. 



flatten the whole system. This tendency of flattening 
the distribution of dwarf galaxies, and claiming the stabil- 
ity of the plane tt, can be considered valid in the limit of 
the linear approximation, i.e. roughly for 150 kpc above 
and below the plane and for a period of time of roughly 9 
Gyr, thus being partially able to imprint the initial proper 
motions of the dwarf galaxies taken into consideration for 
the Local Group. 

Along the directions parallel to the plane, the situation 
is more complicated and described by relations like 



GM\ 



M31 



GMi 



MW 



|^M31 



|»7mw| 



T 12 ri2 + m^l - 2^37)2 + n S T]2 = 

The first two terms in the sum on the left hand side are due 
to the tidal field that has to balance the third term due to 
the centrifugal force, the last term is due to the Coriolis 
effect and the extra term deriving from a non-uniform 
rotation of the system S^. The same kind of relation holds 
for the other coordinates. 

4. Summary, conclusions and consideration on the 
limits of the approach used 

The spatial distribution of the galaxies in the LG is the 
footprint of its formation mechanism, the internal gravi- 
tational interactions among the galaxies, and the gravita- 
tional action of external massive galaxies or galaxy groups 
on the LG members. 

In this paper, we have thoroughly addressed the whole 
subject focusing the attention on the role played by the 
tidal force field exerted by external galaxies or galaxy 
groups on the dwarf galaxies of the LG, excluding those 
that are clearly under the dominant gravitational effects 
of the HGs, in shaping the large-scale distribution of the 
LG galaxies. 

The idea stands on the well known effect of tidal in- 
teractions, which can be expressed as a function of the 
gradient in the gravitational force. While the gravitational 
force never changes sign, its gradient can do so. Moreover, 
while the gravitational force field at any distance from 
the center of mass of a system depends only on the inner 
distribution of matter, the tidal force field does not; it is 
indeed the result of both internal and external distribu- 
tions of matter. The tidal force acting on a body moving 
along a certain direction will pull it away from the ori- 
gin of the reference frame and, at the same time, push it 
along directions perpendicular to the motion toward the 
origin of the reference frame. Therefore a system subjected 
to tidal interactions can undergo not only tidal stripping 
but also tidal compression. In other words, the space dis- 
tribution of galaxies undergoing tidal interactions tends to 
become flat. 

The results of this study can be summarized as follows: 

— The tidal forces can be the cause of the planar distri- 
bution of these dwarf galaxies. We analytically obtain 
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the same numerical results of lPasetto fe Chiosil (|2007l ) . 
In fact, we prove that a planar distribution of all dwarf 
galaxies, excluding those tightly bounded to a HG, is 
compatible with the presence of an external force field. 
T he planar ge o metric al distribution found by 



Pasetto fe Chiosil (|2007f ) was not known to relate 
to the most negative eigenvalue (and associated 
eigenvector) of the tidal tensor. In that sense this 
previous work was partially incomplete. Here we 
have g one deeper into this issue following th e original 
idea of Ravchaudhurv fe Lvnden-Bell (1989) to check 
whether the planar distribution is a mere coincidence 
or the consequence of fundamental laws of mechanics. 
To address this, we first check, using different argu- 
ments, the coincidence between the direction given by 
the vector orthogonal to the geometrical plane and 
that corresponding to the eigenvector with the most 
negative eigenvalue. Second, we analyze the energy 
of the orbital motion of the LG galaxies and find 
that the minimum energy corresponds to a planar 
distribution which is exactly the geometrical plane 7r. 
Therefore the planar distribution is the consequence 
of the long time-scale influence of the tidal forces 
exerted by massive galaxies or galaxy groups external 
to the LG. Finally, we demonstrate that this situation 
has been stable over the past 9 Gyr. 
The equilibrium and stability of the plane is a con- 
sequence of the minimum in the Action and of 
the orbits that c o me fr om this minimum. Although 



Pasetto fe Chiosil (|2007l ) have carefully investigated 



the nature of the Minimum Action (whether local 
or absolute), the uncertainty affecting the orbits de- 
rived from the Action minimization could lead to an 
uncertainty in the energy analysis, the stability of 
the plane and the compression effect in turn here 
evidenced. This problem is still unsolved and can- 
not be resolved at the present time because better 
observational data would be required (proper mo- 
tions, velocities, distance moduli, absolute positions, 
masses, etc.). Hence, the analytical approach pre- 
sented here still requires further support from inde- 
pendent arguments. Along this line of th o ught is the 
study of Ivan der Marel fe Guhathakurtal (|2008h who 
find compatibility betw een their results and what in 
IPasetto fc Chiosil (|2007l ). 

The completeness of the sample of external galaxies 
listed in Table [T] is another factor influencing result s 
developed here as well as in IPasetto fc Chiosil (2007). 
In the previous paper special care was taken to confirm 
the results based on the minimization of the Action 



with an extended catalog from iPeebles et alJ (|200ll ) 



Moreover, we can confirm an excellent concordance be- 
tween the actu al direction of the quadrupole tensor 
eigenvectors in Ravchaudhurv fc Lvnden-Bell (1989) 
and that deri ved here with t he in dependent catalog 
compilation of IPeebles et al.l (|200ll ). Nevertheless the 
continuous discoveries of new LG members suggest 



that the actual census of th e LG galaxies cannot y et 
be considered complete (e.g. Loeb fc Naravanl ( 20081 )). 
The Minimum Action together with the study of the 
equilibrium through the first derivative of the energy 
of the system is a method that can be used to constrain 
the energy of dwarf galaxies with unknown proper mo- 
tions. The missing proper motion prevents us from 
having the complete energy of any given dwarf, but the 
minimization of the action, together with the study of 
the tidal tensor, permits us to handle the derivative of 
the energy in this particular situation. 
We can find another example of a flat structure in 
the super-galactic plane, a slab of roughly = 25 Mpc 
in thickness and of a diameter greater than 110 Mpc 
( Lahav et all 12000). w hich was already prop osed by 
de Vaucouleurs in 1953 (jde Vaucouleurslll953l ). All the 
groups of galaxies used here (Table [I} lie within this 
plane. With respect to this super-galactic plane the 
normal of the plane tt has direction (SGL, SGB) = 
(95°, 69°). This plane has been determined requiring 
that the vector joining the MW and M31 explicitly 
belong to tt. If you attempt to find a best fit determina- 
tion for the plane, say ff, of the overall sample of the LG 
galaxies, without forcing the MW and M31 to belong 
to such a plane, we have, from PC07, (/, b) = (46°, 29°), 
i.e. (SGL, SGB) = (93°, 67°), for the normal to the 
plane ft. Thus, this is even closer to the direction of the 
vector associated with the most negative eigenvalue of 
the tidal tensor but slightly more distant from the di- 
rection of the super-galactic North Pole. The proxim- 
ity of the direction of the normal to tt and the super- 
galactic North Pole cannot yet be claimed as a sig- 
nificant result without further investigation. However, 
it is also not yet possible to study this problem fully, 
given the inc ompleteness of the catalogu e of the nearby 
galaxies (see Karachentsev et all ( 20041 )). The mass es- 
timation together with the distances for the systems 
involved are still the major source of the errors we con- 
sidered. Finally, in none of the external groups taken 
into account in Table [T] can we claim the existence of 
a similar flat distribution of dwarf galaxies that could 
be a first indication of a common external effect act- 
ing on all these groups. It seems that the local dwarf 
galaxies are primarily influenced the external potential 
of nearby groups, not by the larger and more distant 
mass accumulation responsible for the super-galactic 
plane. 

The issue of the spatial distribution of nearby satel- 
lite galaxies bound to their HG remains to be ad- 
dressed. These dwarf galaxies have orbits whose typ- 
ical distance lies inside the dark matter halo of the 
hosting galaxy. Many N-body simulations have shown 
that galaxies in close proximity of a HG are subjected 
to strong collisional interactions. The satellite dwarf 
galaxies may suffer kick-off, interactions and may be 
absorbed into the halo of the HG (MW and M31 
in our case), and may abruptly change the phase- 
space density distribution function. In other words, the 
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large-scale descripti on of the gravitational in teraction 
adopted here and in lPasetto fe Chiosil (|2007l ) does not 
work at the short distance scales of the HG-dwarf satel- 
lite interactions. This will be the subject of a forth- 
coming study, in which simulations of the interaction 
between a HG and its satellites will be investigated in 
the framework of strong collision al dynamics along the 
line of work already init iated by IPasetto et al. ( 2003f ) 



and lPasetto et all (|2009l ) 
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